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ABSTRACT 

In order to gain insight into the physical mechanisms leading to the formation of stars 
and their assembly in galaxies, we compare the predictions of the MOdcl for the Rise 
of GAlaxies aNd Active nuclei (morgana) to the properties of K- and 850/i-selected 
galaxies (such as number counts, redshift distributions and luminosity functions) by 
combining MORGANA with the spectrophotometric model GRASIL. We find that it is 
possible to reproduce the K- and 850^-band datasets at the same time and with a 
standard Salpeter IMF, and ascribe this success to our improved modeling of cooling in 
DM halos. We then predict that massively star-forming discs are common at z ~ 2 and 
dominate the star-formation rate, but most of them merge with other galaxies within 
~100 Myr. Our preferred model produces an overabundance of bright galaxies at 2; < 1; 
this overabundance might be connected to the build-up of the diffuse stellar component 
in galaxy clusters, as suggested by Monaco et al. (2006), but a naive implementation 
of the mechanism suggested in that paper does not produce a sufficient slow-down 
of the evolution of these objects. Moreover, our model over-predicts the number of 
10^° — 10^^ Mq galaxies at z ~ 1; this is a common behavior of theoretical models as 
shown by Fontana et al. (2006). These findings show that, while the overall build-up of 
the stellar mass is correctly reproduced by galaxy formation models, the "downsizing" 
trend of galaxies is not fully reproduced yet. This hints to some missing feedback 
mechanism in order to reproduce at the same time the formation of both the massive 
and the small galaxies. 
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1 INTRODUCTION 

The Lambda Cold Dark Matter (ACDM) cosmology is con- 
sistent with a large body of observations of the large-scale 
Universe (see, e.g., Spergel et al. 2006). Then, the pre- 
dicted hierarchical evolution of Dark Matter (DM) perturba- 
tions subject to gravitational instability provides a standard 
framework to study the formation and evolution of luminous 
structures in the Universe. However, while the cosmological 
framework is fixed with a small uncertainty, several open 
questions, regarding to the formation and evolution of galax- 
ies, arise from the complex evolution of baryons within the 
potential wells of the DM halos. 

Galaxy formation is observationally constrained by 
many multi wavelength surveys of deep fields (e.g. 
C0MB017, Wolf et al., 2001; DEEP2, Davis et al., 2003; 
GOODS, Giavalisco et al. 2004; GEMS, Rix et al. 2004; 
UKIDSS, Lawrence et al., 2006; COSMOS, Scoville et al.. 



2007; etc.). To compare with these datasets, many thou- 
sands of galaxies must be generated by models. This makes 
a straightforward numerical approach problematic (it has 
been attempted, e.g, by Nagamine et al. 2005, Saro et 
al., 2006, Robertson et al. 2007), so that simpler and 
quicker models have been developed, based on sets of recipes 
that address the various processes involved. These "semi- 
analytical" models (see e.g. Somerville et al. 2001, 2004; 
Granato et al. 2004; Menci et al., 2004; Baugh et al. 2005; 
Kang et al., 2005; Bower et al., 2006; Croton et al. 2006; Cat- 
taneo et al., 2006; De Lucia et al., 2006; Monaco, Fontanot 
& Taffoni, 2007) have been tested against an impressive 
number of observational constraints. In this decade-long 
testing process, many specific models have been unsuccess- 
ful in reproducing constraints like the high-mass cutoff of 
the luminosity function (Benson et al. 2003), the level of 
a-enhancement in elliptical galaxies (Thomas 2005; Na- 
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gashima et al., 2005), the redshift distribution of K-band 
sources (Cimatti et al. 2002a) , the surface density of EROs 
(Cimatti et al. 2002b, Daddi et al. 2002; Smith et al. 2001). 

Many of these difficulties have been overcome by later 
versions of the models, but some of them have required 
strong assumptions: for instance, a top-heavy IMF in star- 
bursts is required by Baugh et al. (2005) to reproduce the 
sub-mm counts. In the most recent models (see i.e. Bower 
et al. 2006; Croton et al., 2006) the cutoff of the galaxy lu- 
minosity function and the bimodality of galactic colors is 
obtained only if cooling flows in large halos are quenched 
by AGN feedback, i.e. by the energy emerging from massive 
black holes accreting at a relatively low rate (in the so-called 
radio mode). 

A recently highlighted observational trend of galaxy for- 
mation is the so-called downsizing of galaxies: at variance 
with the hierarchical trend of DM halos, more massive galax- 
ies tend to form their stars earlier and in a shorter period 
than smaller galaxies, which experience more prolonged star- 
formation histories. While this general trend is known to be 
not incompatible with cosmology once stellar and AGN feed- 
back are properly taken into account (see, e.g., Granato et 
al., 2004; Neistcinn, van den Bosch & Dekel, 2006), recent 
observations show that the details, like for instance the (al- 
most) parallel evolution of the star formation density for 
galaxies of diflferent mass (Zheng et al. 2007), are still not 
reproduced. 

Massive elliptical galaxies have always been remarkably 
elusive objects in this regard. The first versions of hierarchi- 
cal galaxy formation models (see, e.g. White 1996) predicted 
that these galaxies form late {z ^ 1) by the merging of al- 
ready assembled discs, while evidence from stellar popula- 
tions (see Matteucci 1996 for a review), the tightness of the 
fundamental plane (Renzini & Ciotti, 1993), the evolution of 
the color-magnitude relation (Kodamaet al., 1998; Blakeslee 
et al, 2003; EUis et al. 2006) and the local M.g2 - cr rela- 
tion (Bernardi et al. 2003) suggested that they formed early 
(z > 2) in a short burst of star formation. Clearly, in a hier- 
archical Universe the age of stars does not need to coincide 
with the assembly age of the galaxy, defined as the time at 
which the most massive progenitor has at least half of the fi- 
nal stellar mass. Massive ellipticals could then be assembled 
late by dry mergers of other ellipticals, so as to preserve the 
oldness of their stellar populations while producing a low 
assembly redshift (De Lucia et al. 2006). This possibility is 
severely constrained by the modest (if any) evolution of the 
high end of the stellar mass function since z = 1 (Cimatti, 
Daddi & Renzini, 2006): while models predict a doubling of 
stellar masses (Dc Lucia & Blaizot, 2007), evidence excludes 
an evolution larger than ~0.2 dcx (as roughly estimated by 
Monaco et al. 2006; see the references therein). Recently, 
Monaco et al. (2006; see also Conroy, Weschler, & Kravtsov, 
2007) have proposed that the formation of a diffuse stellar 
component in galaxy clusters by scattering of stars during 
dry mergers (see, e.g., Murante et al. 2007) may conspire to 
decrease the expected evolution of the high end of the stellar 
mass function. 

The early formation of massive galaxies and their fol- 
lowing (almost passive) evolution arc best constrained by 
deep observations in the sub-mm band, suited to reveal ob- 
scured star formation events at high redshift, and in the 
K-band, suited to probe stellar masses at lower redshift. 



The warm dust present in star forming clouds absorbs 
most of the UV/blue photons emitted by young stars and re- 
processes them to the FIR, becoming the dominant contrib- 
utor in that band. Moreover, the steeply decreasing shape 
of the galactic SEDs from ~100/im to ~1 mm gives a neg- 
ative K-correction that promotes the observation of star- 
bursts in the sub-mm bands up to z ~ 5. Therefore strong 
starbursts at high redshift axe more easily observed in the 
sub-mm than in the optical. The sub-mm emission is mea- 
surable in a few windows, most notably that at 850 /im. Ob- 
servations with the Submillimiter Common-User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Telescope have 
highlighted the presence of a population of high-redshift 
massive starbursts (see, e.g.. Small et al. 1997; Hughes at 
al. 1998; Chapman et al. 2002; Scott et al. 2002) commonly 
interpreted as galaxies forming stars at rates of hundreds if 
not thousands of M© yr~^. Given the poor angular resolu- 
tion of SCUBA images, the identification of optical counter- 
parts is difficult, and can be achieved using interferometric 
images at longer wavelengths. The resulting redshift distri- 
bution is thought to peak at z ~ 2.4 (Chapman et al. 2003, 
2005). 

The if-band is a very good tracer of the stellar mass at 
z < 1.5 (Gavazzi et al. 1996), it is almost unaffected by dust 
extinction, and requires small Jf-corrections that weakly de- 
pend on the morphological type, so it is ideal to follow the as- 
sembly of the bulk of stellar mass. The local Jf-band LF has 
been measured with great accuracy by the 2MASS collabora- 
tion (Cole et al. 2001; Kochanek et al. 2001). In this paper we 
will focus on four K-hand surveys with (photometric or spec- 
troscopic) redshift coverage, (i) The K20 survey (Cimatti et 
al. 2002a, Pozzetti et al. 2003) is a K < 20 (correspond- 
ing to Kab < 21.84) limited sample, covering 52arcmin'^ , 
with a very high redshift completeness (> 90%). (ii) The 
GOODS-MUSIC (Grazian et al. 2006) catalogue is a multi- 
color sample extracted from the deep and wide survey con- 
ducted over the Chandra Deep Field South in the framework 
of the GOODS project. Here we use the Kab < 23.5 limited 
galaxy sample defined in Fontana et al. (2006). This sample 
covers 143.2aT'cmm^; 28% of the galaxies have a spectro- 
scopic redsliifts, used to train photometric redshift estimates 
for all other galaxies, (iii) The UKIDSS Ultra Deep Survey 
(UDS, Dye et al., 2006) galaxy sample is a complete cata- 
log of Kab < 22.5 selected galaxies over O.Gdeg^; for each 
object in the sample a photometric estimate of the redshift 
is provided, (iv) The VIMOS-VLT Deep Survey (VVDS) 
(Le Fevre et al., 2005) is a spectroscopic survey designed to 
measure redshift for ~ 10^ sources selected, nearly randomly 
from a photometric catalogue. Here we consider a /sT-selected 
sample, with either photometric and spectroscopic redshifts 
obtained by Pozzetti et al. (2007) combining a Kab < 22.34 
limited sample (20% redshift completeness) defined over a 
^42arcmin'^ area, with a Kab < 22.84 limited sample (29% 
redshift completeness) defined over a 172arcm.in^ area. 

This paper is the third of a series devoted to describe 
the MOdel for the Rise of GAlaxies aNd Active nuclei (mor- 
gana). With respect to similar models of galaxy formation, 
MORGANA presents a different, more sophisticated treatment 
of the mass and energy flows between galactic phases (cold 
and hot gas, stars) and components (bulge, disc, halo). In 
particular, the process of radiative cooling of the shocked 
gas is treated with a new model (tested against simulations 
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in Viola et al., in preparation), while feedback is inserted 
following the model by Monaco (2004) and galaxy winds 
and super-winds are allowed. The model is described in de- 
tail in Monaco, Fontanot & Taffoni (2007; hereafter paper 
I). The prediction of the properties of the AGN population 
is presented in Fontanot et al. (2006; hereafter paper II), 
and the prediction of the evolution of the stellar mass func- 
tion, together with those of other similar models, has been 
compared to the results inferred from the GOODS-MUSIC 
data in Fontana et al. (2006). As mentioned above, MOR- 
GANA has been used by Monaco et al. (2006) to address the 
lack of evolution of the high end of the stellar mass func- 
tion and its connection with the building of the diffuse stel- 
lar component in galaxy clusters. In this paper we compare 
the predictions of MORGANA to the data mentioned above 
of deep fields in sub-mm (850 fim) and K bands to test to 
what extent the model is able to reproduce the formation 
and assembly of massive galaxies. To this aim, we have com- 
bined MORGANA with the spectrophotometric code GRASIL 
(Silva et al. 1998) that computes the UV to radio SEDs 
of model galaxies, including a three-dimensional bulge-|-disk 
geometry with a two-phase interstellar medium, the radia- 
tive transfer through the dusty ISM, a realistic dust grain 
model, and a direct computation of the dust temperature 
distribution. 

The paper is organized as follows. In section [2Tl1 we de- 
scribe the main properties of the MORGANA model. In section 
I2.2l we describe how we compute luminosity functions, num- 
ber counts and redshift distributions interfacing the output 
of the model with GRASIL. In section[3]we present our results, 
a discussion is given in Section [l] and in Section [5] we give 
our conclusions. Throughout this work we assume, whenever 
necessary, the concordance cosmological model Qa = 0.7, 
Qm = 0.3, Ho = 70 Km s-^ Mpc-^, as = 0.9. 



2 MODEL 

2.1 Galaxy formation model: MORGANA 

MORGANA is described in full detail in Monaco, Fontanot & 
Taffoni (2007), while AGN accretion is described in Fontanot 
et al. (2006). We give here only a brief description, aimed at 
highlighting the main processes included in the model. 



2.1.1 Algorithm 

MORGANA follows the typical scheme of semi-analytic mod- 
els, with some important differences. Each DM halo contains 
one galaxy for each progenitoiQ; the galaxy associated with 
the main progenitor is the central galaxy. Baryons in a DM 
halo are divided into three components, namely a halo, a 
bulge and a disc. Each component contains three phases, 
namely cold gas, hot gas and stars. For each component the 



Each DM halo forms through the merging of many halos of 
smaller mass, called progenitors. At each merging the largest halo 
survives (it retains its identity), the others become substructure 
of the largest one. The main progenitor is the one that survives all 
the mergings. The mass resolution of the box used for computing 
the merger trees sets the smallest progenitor mass, as explained 
in section [2. 1.21 



code follows the evolution of its mass, metal content, ther- 
mal energy of the hot phase and kinetic energy of the cold 
phase. 

The main processes included in the model are the fol- 
lowing. 

(i) The merger trees of DM halos are obtained using 
the PINOCCHIO tool (Monaco et al. 2002; Monaco, Theuns 
& Taffoni 2002; Taffoni, Monaco & Theuns 2002). 

(ii) After a merging of DM halos, dynamical friction, 
tidal stripping and tidal shocks on the satellite (the smaller 
DM halo, with its galaxy at the core) lead to a merger with 
the central galaxy or to tidal destruction as described by 
Taffoni et al. (2003). At each merger a fraction of the satellite 
stars is scattered to the stellar halo component (Murante et 
al. 2007; Monaco et al. 2006). These stars are not associated 
to galaxies but to the intra-cluster light. 

(iii) The evolution of the baryonic components is per- 
formed by numerically integrating a system of equations for 
all the mass, energy and metal flows; this allows not to be 
restricted to linear dynamics. 

(iv) The intergalactic medium infalling on a DM halo 
is shock-heated, as well as the hot halo component of merg- 
ing satellites (which is given to the main halo) and that of 
the main halo in case of major merger (Msat > 0.2 x Mtot). 
Following Wu, Fabian & Nulsen (2001), shock-heating is im- 
plemented by assigning to the infalling gas a specific ther- 
mal energy equal to 1.2 times the specific virial energy, 
— O.Sf/n/MH (where Mh and C/h are the mass and binding 
energy of the DM halo). 

(v) The profile of the hot halo gas is computed at each 
time-step by solving the equation for hydrostatic equilib- 
rium with a polytrophic equation of state and an assumed 
polytrophic index 7p = 1.2. No hot gas is present within a 
cooling radius rcooi , which is set to a vanishingly small value 
at major mergers. 

(vi) The cooling fiow is computed by integrating the 
contribution to radiative cooling of each spherical shell, tak- 
ing into account the heating from (stellar and AGN) feed- 
back from galaxies. Given the importance of cooling for the 
results presented in this paper, we give here some detail of 
the cooling model. If Tgo and pgo are the temperature and 
density of the hot halo gas extrapolated to r = 0, //hotWp 
its mean molecular weight and = rn/cufw the scale radius 
of the halo (of radius rn and concentration Cnfw), then the 
mass cooling fiow Mco.k results: 



Mco.H = 



X X(2/(7p - 1)) 



(1) 



where the integral X(a) is defined as J^'^"'"^^ {1 — a[l — ln(l + 

t)/t]}'^t^dt, with a = [3Tvir(7p-l)Cnfw(T + C„fw)]/{7pr90[(l + 

Cnfw) ln(l + Cnfw) — Cnfw]} (T^vir bciug the vlrial temperature 
of the halo). The cooling time tcooi.o is computed using the 
central density (the density gradient is taken into account 
by the integral T) and the temperature at rcooi (the temper- 
ature gradient is neglected): 

_ 3fcrg (rcooi )photn^p 

Icool.O n„ /A r \ ^ ' 

^pgO{l\-cool — t heat J 

Here Acooi is the metal-dependent Sutherland & Dopita 
(1993) cooling function and the heating term Fheat is com- 
puted assuming that the energy fiow -Bhw.H fed back by the 
galaxy is given to the cooling shell: 
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47rrfX(2/(7p - 1)) 1 p<,o 



(3) 



Whenever Fhcat > Acooi the coohng flow is quenched. The 
coohng radius rcooi is treated as a dynamical variable whose 
evolution takes into account the hot gas injected by the cen- 
tral galaxy (Afhw.n): 



A/c. 



^cool = -, . (4) 

This equation is valid if pressure is balanced at rcooi, an 
assumption which clearly does not hold in general. We mimic 
the pressure force acting on mass shell at rcooi as follows: 

^cool ~ rcooi ~ Cs (5) 

where Cs is the sound speed computed at rcooi. This recipe is 
different from that used in most other semi-analytical mod- 
els, where a cooling radius is computed by inverting the 
cooling time as a function of radius. Viola et al. (2007) have 
compared analytic cooling models to N-body -\- hydro sim- 
ulations showing that, while the "classical" cooling model 
significantly underestimates the amount of cooled mass, the 
present cooling model gives a very good fit. 

(vii) When the hot halo phase is heated by feedback 
beyond the virial temperature, it can leave the DM halo 
in a galactic super-wind. A similar thing happens to the 
cold halo gas when it is accelerated by stellar feedback. To 
compute the time at which the ejected gas falls back into a 
DM halo, its merger history is scrolled forward in time until 
the circular velocity is larger than the (sound or kinetic) 
velocity of the gas at the ejection time. 

(viii) The cooling gas is let infall on the central galaxy 
on a dynamical time-scale (computed at rcooi)- It is divided 
between disc and bulge according to the fraction of the disc 
that lies within the half-mass radius of the bulge. 

(ix) The gas infalling on the disc keeps its angular mo- 
mentum; disc sizes are computed with an extension of the 
Mo, Mao & White (1998) model that includes the contribu- 
tion of the bulge to the disc rotation curve. 

(x) Disc instabilities and major mergers of galaxies lead 
to the formation of bulges. We also take into account a pos- 
sible disc instability driven by feedback. In minor mergers 
the satellite mass is given to the bulge component of the 
larger galaxy. 

(xi) Star formation and feedback in bulges and discs are 
inserted following the model of Monaco (2004). According 
to that model, the regime of stellar feedback in a galaxy 
depends mainly on the density and vertical scale-length of 
the galactic system. In thin systems, like spiral discs with 
gas surface density Ecom.d and fraction of cold gas /coW.d, 
the timescale for star formation tt.D is predicted to be: 

\ -0.73 0.45 
i^cold.D \ / /cold.D 



U.u = 9.1 



1 Mq pc- 



/ Jcold.D \ 
V 0.1 ) 



Gyr 



(6) 



Due to the correlation of /coM.d and Scold, d (galaxies with 
higher gas surface density consume more gas), this relation 
is compatible with the Schmidt law. Thick systems like star- 
forming bulges (or mergers) are dominated by transients 
which are very difficult to model, so the straightforward 
Schmidt law is used: 



J cold, B 



1 M 



PC- 



Gyr 



(7) 



In both cases, hot gas is ejected to the halo (in a hot galactic 
wind) at a rate equal to the star-formation rate (as predicted 
by Monaco 2004), though massive bulges with circular ve- 
locity Vb > 300 km s~^ are able to bind their hot phase 
component. The thermal energy of this re-heated gas is not 
scaled to the DM halo circular velocity but to the energy of 
exploding SNe. 

(xii) In star-forming bulges cold gas is ejected in a cold 
galactic wind by kinetic feedback due to the predicted high 
level of turbulence driven by SNe (this process is described 
in full detail in paper II, Section 2.2). 

(xiii) Accretion of gas onto massive black holes (start- 
ing from small seeds present in all galaxies) is connected to 
the ability of cold gas to loose angular momentum by some 
(unspecified) mechanism driven by star formation. This is 
explained in full detail in paper II. 

(xiv) Because accretion onto black holes is triggered 
by star formation, AGN feedback can quench cooling flows 
only after their start. Alternatively, the cooling flow can be 
quenched when a fiducial energy criterion is met, without 
letting any gas fall to the galaxy and form stars. This "forced 
quenching" procedure has been used in paper II and will be 
discussed in a forthcoming paper. 

(xv) Metal enrichment is self-consistently modeled in 
the instantaneous recycling approximation. 

2.1.2 Runs 

All models are based on the same pinocchio run we intro- 
duced in paper I, a 512'^ realization of a ISOAfpc comoving 
box {h = 0.7). The mass particles is 1.0 x 10^ M0, and the 
smallest halo we consider contains 50 particles, for a mass 
of 5.1 X 10^" M0. The branches of the DM halo merger 
trees are tracked starting from a mass of 10 particles, corre- 
sponding to 1.0 X 10^'^ Mq; this is the mass of the smallest 
progenitors. For sake of comparison, the corresponding val- 
ues for the Millennium Simulation (Springel et al. 2005) are 
1.2 X 10^ for the particle mass and 2.5 x 10^° for the smallest 
progenitor. We have tested the overall stability of our results 
by running the model over two other 512'^ boxes of size 200 
and 100 Mpc (see paper I for more details). 

The stellar mass of the typical galaxy contained in the 
smallest DM halo at 2: = is 3 x 10* M0 . In order to estimate 
the completeness limit for the stellar mass function we con- 
sider an higher resolution pinocchio run (512^ realization 
of a IQQMpc comoving box, see also paper I, appendix B, for 
a discussion about the numerical stability of the model). We 
study the typical stellar mass of the galaxies belonging to 
2.5 - 5.1 X 10^° M0 DM halos and we estimate a complete- 
ness limit of 4.5 x 10* Mq. For each run we compute the 
evolution of (up to) 300 trees (i.e. DM halos at 2; = 0) per 
logarithmic bin of halo mass of width 0.5 dex. This implies 
that while all the most massive halos are considered, smaller 
halos are randomly sparse-sampled. To properly reconstruct 
the statistical properties of galaxies we assign to each tree a 
weight Wtreo equal to the inverse of the fraction of selected 
DM halos in the mass bin. 

A standard parameter choice was presented in paper 
I. However, it was noticed there that spiral discs tend to 
be more compact than in the real Universe. A Schmidt- 
Kennicutt-like law (equation |6]) is used to compute star 
formation rates, so this results in higher surface densities. 
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shorter star formation timescales and lower gas fractions at 
z = 0. Forcing DM halo concentrations to lower values allows 
to alleviate this problem, while reproducing the zero-point 
of the TuUy-Fisher relation and preserving a good fit of the 
Schmidt-Kennicutt law. We then scale concentrations so as 
to take a given value for a lO^'^ Mq DM halo at z = 0. 
The value was proposed to be 4 in paper I; however, we 
noticed that this choice tends to lower the stellar mass func- 
tion at the knee, so we use the slightly higher value of 6 
as a good compromise. Moreover, following Monaco et al. 
(2006) we scatter to the halo stellar component 40% of the 
stellar mass of satellites at each galaxy merger. Finally, we 
implement AGN feedback by adopting the "forced quench- 
ing" procedure of paper II; this choice has little influence on 
the results presented here, and will be discussed in a forth- 
coming paper. Also, we have checked that the inclusion of 
quasar-triggered galaxy winds, necessary to reproduce the 
accretion history of massive black holes (paper II), does not 
change the qualitative behavior of the models. 

These changes of parameters do not influence in any 
way the conclusions drawn in this paper. In fact, we are not 
proposing a "best fit" model of the galaxy population; on 
the contrary, we will show that MORGANA is able to repro- 
duce many observables, but its agreement with data breaks 
as soon as one tries to reproduce some relevant aspects of 
the "downsizing" trend, so that a global flne-tuning of the 
model is not possible. The most important point of such 
an investigation is to understand the limits of these models 
and the origin of the disagreement with data. Fine-tuning to 
fit specific datasets, like the local luminosity functions, will 
be required in other contexts, for instance to create galaxy 
mock catalogues. 

2.2 SED model: GRASIL 

For each galaxy modeled by MORGANA we compute the cor- 
responding UV to radio SED with GRASIL. The details of the 
code are given in Silva et al. (1998) (and the subsequent up- 
dates and improvements in Silva 1999; Granato et al. 2000; 
Bressan, Silva, Granato 2002; Panuzzo et al. 2003; Vegaet al. 
2005), while we summarize here the main features, (i) Stars 
and dust are distributed in a bulge (King profile) -I- disk (ra- 
dial and vertical exponential profiles) axisymmetric geome- 
try, (ii) We consider the clumping of both (young) stars and 
dust, through a two- phase interstellar medium with dense 
giant molecular clouds (MCs) embedded in a diffuse ("cir- 
rus") phase, (iii) The stars are assumed to be born within 
the optically thick MCs and gradually to escape from them 
as they get older on a time-scale tcsc- This gives rise to an 
age-selective extinction, with the youngest and most lumi- 
nous stars suffering larger dust extinction than older ones, 
(iv) The dust consists of graphite and silicate grains with a 
distribution of grain sizes, and Polycyclic Aromatic Hydro- 
carbons (PAH) molecules. In each point of the galaxy and 
for each grain type the appropriate temperature is computed 
(either equilibrium T for big grains or probability distribu- 
tion of temperature for small grains and PAHs). The de- 
tailed PAH emission spectrum has been updated in Vega et 
al. (2005) based on the Li & Draine (2001) model, (v) The 
radiative transfer of starlight through the dust distribution 
is computed yielding the emerging SED. The simple stellar 
population (SSP) library (Bressan, Granato, & Silva 1998; 



Table 1. Adopted values for the GRASIL parameters that are not 
provided by MORGANA. 



Parameters Values 

tesc 10'^ yr 
/mc 0.5 
MMc/rlic M0/(16 pc)2 
f^^K 0-1 
m 0-1 

Bressan et al. 2002) includes the effect of the dusty envelopes 
around AGB stars, and the radio emission from synchrotron 
radiation and from ionized gas in HII regions. 

The GRASIL parameters that are not provided by MOR- 
GANA are summarized in Table [T] and are described in the 
following, (i) The escape time-scale of young stars for the 
parent MCs is set to tosc = 10^ yr. This is intermediate 
between the values found by Silva et al. (1998) to well de- 
scribe the SED of spirals (~ a few Myr) and starbursts ('^ 
a few 10 Myr), and is of the order of the estimated de- 
struction time scale of MCs by massive stars (e.g. Monaco 
et al 2004b). (ii) The gas mass predicted by MORGANA is 
subdivided between the dense and diffuse phases, fixing the 
fraction of gas in molecular clouds /mc to 0.5. The resulting 
SEDs are not much sensitive to this choice, (iii) The mass 
of dust is obtained by the gas mass and the dust to gas 
mass ratio S which is set to evolve linearly with the metal- 
licity {6 = 0.45 Z). (iv) The optical depth of MCs, driving 
their spectrum, depends on the mass and radius set for MCs 
through T cc S Mmc/''mci with Mmc = 10® -^0) ''mc = 16 
pc; (iv) the bulge and disk scale radii for stars and gas are 
given by MORGANA, while the disk scale heights, h"^ and h'^ 
for stars and dust respectively, are set to 0.1 the correspond- 
ing scale radii. In addition, the dust grain size distribution 
and abundances are set to match the mean Galactic extinc- 
tion curve and emissivity (as in Silva et al. 1998 and Vega 
et al. 2005) and are not varied here. 

2.3 Interfacing MORGANA with GRASIL 

The output of MORGANA consists, for each galaxy, in a time 
sampling of the main dynamical variables of the model: for 
each component (halo, bulge, disc) the code issues mass, 
kinetic energy and metal mass of cold gas; mass, thermal 
energy and metal mass of hot gas; mass and metals of 
stars; average (over the time bin) and punctual (at the 
end of the time bin) star formation rate (these are given 
only for bulge and disc). For each galaxy it gives also 
expelled mass and metals, black hole mass, cooling radius, 
bulge and disc radii and velocities, punctual values of the 
accretion rate onto the black hole. Star formation histories 
are reconstructed following all the exchanges of stellar 
matter between galaxies and galaxy components, so that 
they refer to the stars contained in the component at a 
given time and not to the stars formed in that component. 

Of this information GRASIL uses the star formation and 
cold gas metallicity history for bulge and disk, the mass of 
gas in the two components at the required galactic age, and 
the stars and gas scale radii at the same age. The time sam- 
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Figure 1. X-band Luminosity Function. Data refer to the ob- 
servations of Cole et al. (2001) and Kochanek et al. (2001). Solid 
line refers to MORGANA prediction. 

pling of these quantities is 0.1 Gyr. GRASIL resamples the 
star formation histories on a much finer time grid to opti- 
maUy account for the short lifetimes of massive stars. How- 
ever, the width of the time bin of MORGANA is rather large, 
so the number of massive stars can be rather inaccurate. We 
then split the last time bin into a sub-bin of 10 Myr, to which 
we assign the punctual value of the star-formation rate at 
the end of the bin, and a larger, earlier one of 90 MyiQ, to 
which we assign a star formation rate such that the inte- 
gral in the two sub-bins gives the correct final amount of 
stars. Using test star-formation histories, we have verified 
that this sampling allows to reproduce the magnitudes to 
within <0.1 mag; indeed, even the strongest star-formation 
events do not have associated time-scales much smaller than 
the time bin. 

It is worth mentioning that running the spectro- 
photometric code on model galaxies is the main bottleneck 
of the computation; it is then necessary to device strategies 
to optimize the computations by estimating the minimal 
number of galaxies needed to have a reliable result. 

2.4 Luminosity functions, number counts and 
redshift distributions 

To simulate a deep field, the computation of galaxy SEDs is 
performed on the box at several output times. As explained 
in paper I, the sparse-sampling of trees (section 12.1.2^ re- 
sults in an over-sampling of small satellites with respect to 

^ As the age of Universe in the assumed cosmology is 13.47 Gyrs, 
the time bin corresponding to z = is smaller than 0.1 Gyr; in 
this case we make the earlier bin consistently smaller than 90 
Myr. 



central galaxies of similar stellar mass. We correct for this 
over-sampling by further sparse-sampling the satellites as 
follows. First, we construct from the results of the run at 
a = an average curve of mass of the central galaxy as 
a function of DM halo mass; second, we randomly sparse- 
sample the satellites with a probability equal to the ratio 
between the weights lutrco of the tree the satellite belongs 
to and that of the DM halo whose central galaxy has on 
average the same mass as the satellite. The inverse of this 
probability is a new weight, Wsat- Central galaxies are all se- 
lected and assigned a unity weight. If a galaxy is destroyed 
by mergers (or tides) then its stellar mass at the destruction 
time is used to compute Wsat- This procedure gives a roughly 
constant number of galaxies in logarithmic intervals of mass. 
Because each galaxy is present in many time bins, the num- 
ber of selected galaxies is typically very high. To limit this 
number we introduce a third weight uigai as the inverse of 
a further sampling factor, equal for all galaxies. The first 
two samplings (of merger trees and galaxies) are both com- 
puted at 2 = 0, but the weights wtrcc and Wsax assigned 
to the galaxies are used at any redshift. This is done with 
no loss of generality, as a fair reconstruction of luminosity 
functions or number counts only requires that the weights 
are used consistently with the sparse-sampling procedure; in 
other words, we only need to require that the properties of 
a sparsely sampled population are weighted by the inverse 
of the sampling probability, whatever the redshift at which 
the sampling is performed. 

In order to ensure a smooth description of the redshift 
evolution of the properties of our simulated galaxies, we use 
the following procedure. At every integration time-bin, cor- 
responding to a redshift z(t), at the end of integration we 
apply the three sampling procedures to all the galaxies in 
the box and select a subsample. It is worth to notice that 
the random sparse sampling ensures that a different subsam- 
ple of galaxies is considered at each z{t), so that galaxies at 
different redshifts are not always the same galaxies seen at 
different times. We define a fourth weight Wtimc as the ratio 
between the width of the time sampling (O.lGj/r) and the 
time interval span by l^ox at z{t). We also compute the an- 
gle subtended by a square of comoving side Ihox at z{t). We 
then compute the SEDs of model galaxies using GRASIL and 
we use them to estimate absolute and apparent magnitudes 
(both in the Vega and AB system) and infrared fluxes. We 
collect these informations in catalogues. 

Luminosity functions, number counts, and redshift dis- 
tributions of magnitude-limited samples are then computed 
by performing weighted sums over the galaxies using the 
product of the four weights defined above. Luminosity func- 
tions at high redshift are computed over the same redshift 
intervals as in the observations, while number counts and 
redshift distributions of magnitude- limited samples are com- 
puted using galaxies starting form z ~ 6. 



3 RESULTS 
3.1 iC-band 

The resulting K-hand LF at z = (in the Vega sys- 
tem) is compared in fig. [1] to that obtained using the local 
2MASS sample (Jarrett et al. 2000), limited to Kvega < 13.5 
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Figure 2. Source Number Counts in the ii'-band. Data refer to 
observations as listed in the figure. SoUd Une refers to MORGANA 
prediction. 
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Figure 3. Redshift source distribution compared to K20 sample 
(histogram, Cimatti et al. 2002a). Solid line refers to MORGANA 
prediction. 



(Kochanek et al. 2001), and the combination of the 2MASS 
and 2dF samples (Cole et al. 2001). Being the K-band lu- 
minosity a good tracer of stellar mass, this comparison is 
analogous that of figure 7 of paper I, where the stellar mass 
function was compared to that inferred by data (obtained 
also with the same 2MASS+2dF sample). In agreement with 
that result, the z = K-hand LF shows an overestimate of 
both the bright (massive) end and the faint (low-mass) end 
of the LF. In particular, the overestimate of the high tail is 
much more relevant here, with respect to the estimated mass 
function. We find that this is in part due to a discrepancy in 
the adopted Mt/Lx ratios when passing from the luminos- 
ity function to the mass function. The average M^/Lk ratio 
of the MORGANA-I-GRASIL model galaxies is ~ 0.7 Mq/Lx©, 
while for instance Cole et al. (2001) adopt an higher average 
value, M^/Lk ~ 1.32 Mq/Lkq, typical of very old stellar 
populations (see e.g. figure 24 in Maraston 2005). In spite of 
these discrepancies, the model is able to reproduce correctly 
the overall normalization of the LF. 

Fig. [2]shows the predicted K-hand number counts com- 
pared to data available in the literature (Mobasher et al. 
1986; Glazebrook et al. 1993; Gardner et al. 1993, 1996; 
Djorgovski et al. 1995; Moustakas et al. 1997; Szokoly et al. 
1998; Kochanek et al. 2001; Maihara et al. 2001; lovino et al. 
2005; Kong et al. 2006; Grazian et al. 2006). The models fit 
well the data in the range 15 < Kvega ^ 22. The excess at 
bright fluxes is a signature of the overestimate of the LF at 
the bright end, while the excess at faint fluxes is commented 
below. 

Number counts do not strongly constrain the model 
unless redshift distributions of magnitude-limited samples 
are available. We then compare our results with redshift 
distributions obtained from the K20 (fig. |3l Vega system), 



GOODS-MUSIC (fig.H AB magnitudes) and VVDS (fig.El 
AB magnitudes) catalogues. For sake of clarity we show 
both the differential and the cumulative distributions. We 
compute the error on the cumulative distributions using a 
bootstrap technique based on 1000 mock redshift catalogues 
drawn using the observed redshift distribution. The agree- 
ment of MORGANA prediction with K20 observations (rela- 
tive to Kab < 21.84) is good: we predict a total number 
of 435 objects against 480 ± 22 observed. The position of 
the peak of the distribution and the long tail of galaxies at 
z > 1.5 are both recovered; we relate the disagreement in the 
cumulative distribution at low redshift to the known cluster 
at 2 ~ 0.7; some excess at z ~ 0.5 can be again connected 
to the excess of bright galaxies in the K-band LF. The com- 
parison with the deeper but less wide GOODS-MUSIC sam- 
ple (Fig. |4]) confirms the good agreement at brighter limit 
magnitudes (301 predicted vs 262 observed galaxies with 
Kab < 20; 722 predicted vs 623 observed galaxies with 
Kab < 21). At fainter limit magnitudes the comparison 
highlights an excess of model galaxies (1687 predicted vs 
1258 observed galaxies with Kab < 22; 3572 predicted vs 
2603 observed galaxies with Kab < 23). This excess is due 
not only to the small local excess noticeable in the z = 
LF (Fig. [1} , but to a generalized excess of faint sources at 
all redshifts. Fontana et al. (2006) showed that the over- 
prediction of small galaxies (Af ~ 10^" Mq) at z ~ 1 is a 
common problem of galaxy formation models; here we see 
that the problem is probably present since high redshift. Fi- 
nally, we also compare MORGANA predictions with the larger 
VVDS sample (Fig. [5] we combined model predictions for 
the two sub-samples separately in the same way as obser- 
vational data): we obtain again a good agreement with ob- 
servations (9822 predicted objects against 10160 observed); 
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Figure 4. Redshift source distribution compared to GOODS-MUSIC sample (histograms, Fontana et al. 2006). Left Panel: redshift 
distributions. Right Panel: cumulative distributions. Solid line refers to MORGANA predictions. 
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Figure 5. Redshift source distribution compared to VVDS sam- 
ple (histogram, Pozzetti et al. 2007). Solid line refers to MORGANA 
prediction. 

except for a mild underestimate, mostly due to the smaller 
number of 2 ~ 1 — 2 objects with respect to observations. 

A clearer view is obtained by considering the evolution 
of the K-band luminosity function with redshift (Fig. 
We compare our model with the K20 (Pozzetti et al. 2003; 
Vega system) and UDS (Cirasuolo et al. 2006; AB mag- 



nitudes) data; the K20 estimate is based on spectroscopic 
redshifts but the sampled area is small, the UDS sample 
covers a larger area of the sky but is based on photometric 
redshifts. The overall agreement between MORGANA and the 
two datasets is very good. We can notice some overestimate 
of the faint-end at z ~ 0.5, which is not as visible as in 
Fontana et al. (2006). Also, the lower redshift bins of the 
wider UDS sample show that the excess at the bright end 
builds up at 2 < 1. 

From the analysis of fig. [1] to fig. |6] we conclude that 
MORGANA is able to reproduce the assembly of the bulk of 
the stellar mass, which is mostly in place already at 2: ~ 1 — 
2; however, the biggest model galaxies continue growing in 
mass at 2 < 1, though this growth is partially compensated 
for by the loss of stars to the halo by gravitational scattering. 
Moreover, the progressive build-up of smaller objects is not 
correctly reproduced; small objects tend to be too many at 
2 ~ 1 (and too old at 2 ~ 0) in the model. These findings 
are in agreement with the analysis presented in Fontana et 
al. (2006). 

3.2 850^m band 

We now test whether the stars visible at 2 < 2 in the K-band 
were formed in much smaller chunks that were later assem- 
bled together, or in big massive starbursts at 2 > 2. This 
is best tested by comparing the model to number counts at 
850/im, where massive starbursts at 1 < 2 < 5 are easily 
observed. Fig. [7] shows our prediction compared to SCUBA 
observations (Hughes et al. 1998; Barger et al. 1999; Eales 
et al. 2000; Webb et al. 2002; Borys et al. 2002, 2003; Chap- 
man et al. 2002; Cowie et al. 2002; Scott et al. 2002; Small 
et al. 2002; Coppin et al 2006; Scott et al 2006). In order to 
avoid contamination due to low-redshift sources, we consider 
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Figure 6. LF Rcdshift Evolution. Left panels: K20 sample, data 
points taken from Pozzetti et al. (2003). Mid and Right Panel: 
UDS sample, data points taken from Cirasuolo et al. (2006). Solid 
line refers to MORGANA predictions. 
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Figure 7. SSO^im Source Number Counts. Data refer to observa- 
tions as listed in the figure. Solid line refers to MORGANA predic- 
tion. 



Figure 8. Predicted Redshift Distribution for 850/im sources 
with S > 0.2mJy. 

here a subsample with z > 0.5. We stress that this predic- 
tion is done assuming a standard Salpeter IMF. Fig.[S]giv6s 
the redshift distribution of the objects with flux > 0.2 mJy, 
which corresponds to the faintest point in fig. [T] It is peaked 
at z ~ 2.0, in qualitatively agreement with the observa- 
tions of Chapman et al. (2003,2005) of the brightest SCUBA 
sources (redshift distribution peaked at z ~ 2.4 with a quar- 
tile range of 1.9 < z < 2.8). The model reproduces well the 
data, with a possible modest overestimate (underestimate) 
at faint (bright) fluxes. However, a deeper analysis of the 
predictions shows that the population of brightest objects 
(> 5mJy) at z > 2 (Chapman et al., 2005; Aretxaga et al., 
2007), is missing: we are able to reproduce the bulk of the 
heavily star-forming population at 2 ~ 2, but not the most 
luminous objects. 



4 DISCUSSION 

Sub-mm counts have longly been an elusive piece of evi- 
dence to fit for semi-analytic models of galaxy formation 
(e.g. Guiderdoni et al. 1998; Devriendt & Guiderdoni 2000; 
Granato et al. 2000; but see Baugh et al. 2005). Conversely, 
in the long tuning and debugging process of the MORGANA 
code, we have never found any difficulty in fitting the bulk 
of number counts using a conservative choice for the IMF. 
The difference can be ascribed to the different cooling model 
that MORGANA is implementing. In a forthcoming paper (Vi- 
ola et al. 2007) we compare analytic cooling models to N- 
body-fhydro simulations in the simplest case of an isolated 
halo with its hot gas component initially in hydrostatic equi- 
librium, and no feedback from star formation or AGN. The 
classical model of cooling, based on the computation of a 
cooling radius (White & Frank 1991) is found to underesti- 
mate the amount of cooled gas, especially in the first stages 
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Figure 9. Left panel: star formation rate (bulgc+disc) versus cooling rate for central galaxies at z = 2.5. Right panel: ratio of the cooling 
flows as predicted by the classical model and as used by MORGANA versus total star formation rate for the same galaxies as above. In 
the right panel many points, especially at low star-formation, lie out of the plot because the classical model predicts that all the gas has 
cooled, while MORGANA retains some gas through feedback and cosmological infall. 



of cooling, while the model implemented in MORGANA and 
briefly described in Section 12. f I gives a very good fit. All 
the analytic models of cooling tend to give similar results 
at later times. The use of a better cooling model acceler- 
ates the accumulation of cold gas in the halos, giving rise 
to stronger starbursts. The typical mass of the halos that 
contain the biggest starbursts is high enough to ensure the 
formation of a hot halo component through shocks (see, e.g., 
Keres et al. 2005), so that the results based on hot gas halos 
in hydrostatic equilibrium apply to this case. 

To support this interpretation, we run again the model, 
computing at each halo major merger the gas profile (as de- 
scribed in Section [2.11) and from it the cooling time tcooi{r) 
as a function of radius. The classical cooling radius rc{t) 
is the inverse of the function tcooi(''). The cooling flow of 
the classical cooling model is then computed as Mdassicai ~ 
A-rvrQpg (rc)rfrc / dt. It is not possible to use the classical cool- 
ing flow directly in the model, because the cooling model in- 
cludes the injection of energy by feedback, which would then 
be absent. Besides, a generalization of the model to include 
classical cooling would require deep changes in the code and 
a re-calibration of the parameters to reproduce local observ- 
ables, which is well beyond the interest of the paper. We then 
simply compare the classical cooling flow with that actually 
used in the MORGANA model. Firstly, in figure (9] left panel, 
we show MORGANA cooling fiow versus star formation for 
all the central galaxies present in the box at z = 2.5. The 
stronger starbursts are associated with the stronger cooling 
fiowqj. Second, we show for the same halos (right panel) the 

^ A few massive starburst are associated with no cooling flow and 
then lie outside this relation; these are mergers where feedback has 



ratio of classical and MORGANA cooling fiows versus star for- 
mation. For massive starbursts (A/^f > 100 M0 yr~^) the 
classical model predicts cooling fiows a factor of 5-10 lower, 
which would presumably result in correspondingly lower star 
formation rates. The relation changes at smaller star forma- 
tion rates, where feedback has been able to suppress the 
MORGANA cooling fiow in many casefl The success of MOR- 
GANA should then be ascribed to the different modeling of 
cooling. 

With our model we can also investigate the physical 
conditions of SCUBA sources. In fig. llOb , we show the star- 
formation function of galaxies a,t z — 2.5, divided into bulges 
(mergers) and discs. In this case, we use the punctual value 
of the star formation rate computed at the end of the time- 
bin. Interestingly, mergers dominate only at the highest star 
formation rates, > 200 M© yr~^. This implies that most of 
our starbursts are triggered by cooling/infall more than by 
mergers, in line with the findings of violently star-forming 
discs at 2: ~ 2 (Genzel et al. 2006). However, this does 
not imply that massive starbursts are associated with spiral 
galaxies. In Fig. llOb the star formation rate is computed 
as the total amount of stars formed in the 0.1 Gyr time 
bin and eventually found (at the end of the time bin) in a 
bulge or disc component, divided by the width of the time 
bin. This time bulges clearly dominate the star-formation 
function. This shows that such massively star-forming discs 



already quenched cooling, or where cooling has already deposited 
most mass in the galaxy 

* Many points at low star-formation lie out of this plot because 
the classical model predicts no cooling, while in MORGANA gas is 
still present due to feedback and cosmological infall. 
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Figure 10. Star-formation rate function at 2 = 2.5 for discs (blue dotted line), bulges/mergers (red dashed line) and total (black 
continuous line). The star formation in panel (a) is the punctual value at the end of the integration over the time bin, in (b) is the 
average in the bin (of width 0.1 Gyr) relative to the stars contained in bulges and discs at the end of the integration. 



merge into bulges (disc instabilities are not relevant at this 
redshift) in less than 0.1 Gyr. Such short time-scales guar- 
antee high levels of a-enhancement, so the stellar popula- 
tions formed in these starbursts will be typical of early-type 
galaxies. We conclude that the cooling/infall domination of 
high-redshift starbursts does not change the conclusion that 
these events are responsible for the formation of the stars 
found today in elliptical galaxies. 

The overall agreement between model and data is in 
line with recent results by De Lucia et al. (2006), Bower 
et al. (2006), Croton et al. (2006) and Kitsbichler & White 
(2006). However, there are three main points of disagree- 
ment between our model and the data. First, the evolution 
of the most massive galaxies at 2: < 1 is too strong, resulting 
in an excess of bright galaxies in the K-band. It is easy to 
absorb this excess at z = by tuning the parameters, but 
this would be done at the expense of a poor fit of the stellar 
mass function and if-band LF at z ~ 1. As demonstrated 
by Monaco et al. (2006), this evolution is driven by galaxy 
mergers, so no feedback recipe can solve it. Our implemen- 
tation of the suggestion of Monaco et al. (2006) to scatter 
40% stars to the diffuse component of galaxy clusters at 
each merger goes in the right direction but is largely insuf- 
ficient to suppress this evolution. A much higher scattered 
fraction, like ~80% of the stellar mass of the satellite, would 
give a better results (as argued also by Conroy et al. 2007 
and Renzini 2007). However, such an extreme value applied 
to all mergers and all redshift does not allow to reproduce 
the ~ 1% fraction of diffuse stars in the Milky Way together 
with the ~ 10% fraction in Virgo and the ~ 40% fraction 
in massive clusters. We conclude that, while the mechanism 
is promising and deserves further attention, an easy and 
straightforward implementation does not work; this mecha- 



nism should be selectively efficient in the most massive DM 
halos at low redshift. 

Second, the model is not able to reproduce the most 
luminous SCUBA galaxies at redshift 2: ~ 2, which corre- 
spond to the strongest starbursts. This also hints to an in- 
sufficient "downsizing" in our predicted galaxy population. 
In the most massive galaxies the fraction formed at high red- 
shift in episode of intense star formation is still too small, 
and the fraction of stars accreted at later times is still too 
large, with respect to the observed trends. Given the success 
in reproducing the bulk of the SCUBA population, a detailed 
investigation on the predictions of the cooling/infall model 
in the more massive halos is promising. Also the choice of 
a less conservative IMF, i.e. the Kroupa IMF, can help to 
solve the problem, thanks to the larger number of massive 
stars predicted at each stellar generation. 

Third, we confirm the finding of Fontana et al. (2006) 
that our model produces too many faint objects at z ~ 1; in 
other words, our smaller galaxies are too old. Coupled with 
the excessive evolution of the most massive galaxies, this 
finding shows that, while reproducing the average build-up 
of galaxies, the observed trend of less massive galaxies being 
on average younger is not properly reproduced by MORGANA. 
Fontana et al. (2006) have shown very convincingly that this 
problem is shared by many galaxy formation codes, either 
numerical or semi-analytic, so our results can be considered 
as the typical outcome of the ACDM cosmogony, given the 
present understanding of the physics of galaxy formation. 
Figure [TT1 shows a comparison of the evolution of the stellar 
mass function as inferred by the GOODS-MUSIC sample 
and as predicted by morgana. Clearly, the "downsizing" 
trend of more massive galaxies evolving very slowly a.t z < 1, 
while the population of less massive galaxies builds up, is 
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Figure 11. Evolution of the stellar mass function according to MORGANA (panel a) and GOODS-MUSIC (panel b). The latter mass 
functions are drawn roughly in the same mass interval where data are available. 



not reproduced. While the evolution of the high-mass end is 
driven by mergers, a delay in the build-up of faint galaxies 
should be caused by some source of feedback. 



5 CONCLUSIONS 

This paper is the third of a series devoted to presenting 
MORGANA. We have demonstrated that the model is able to 
follow the build-up of the bulk of stars, more precisely to 
reproduce the early assembly and late, almost-passive evo- 
lution of massive galaxies. To this aim we have combined 
MORGANA with GRASIL and compared predictions with ob- 
servation in the submm (the 850/im channel), especially sen- 
sitive to the strongest and most obscured episodes of star 
formation, and in the NIR (especially in the if-band), most 
sensitive to the stellar mass. Overall consistency between 
model and observations has been obtained using SCUBA 
counts at 850/im, number counts in the if-band, redshift 
distribution of A"-limited galaxy samples, and redshift evo- 
lution of the A'-band LF. The importance of this result is 
strengthened by the use of a standard Salpeter IMF (e.g. a 
very conservative choice) along the whole redshift interval. 
We ascribe this agreement mostly to our improved model 
for cooling/infall, that correctly reproduces the results con- 
trolled numerical experiments. 

We predict that most star formation at high redshift is 
not stimulated by starbursts but is due to the strong cool- 
ing/infall flows that take place at early times. This does 
not imply that such discs are close analogues of local spi- 
ral galaxies, as their gas surface density is very high, their 
star formation rate is typical of starburst galaxies and most 
of them are expected to merge within less than 0.1 Gyr. 
We then predict that gas-rich discs, characterized by star- 



formation rates up to 100 — 200 M© yr , should be very 
abundant at z ~ 2, in line with the observations of Genzel 
et al. (2006). 

Despite these successes, our model does not reproduce 
the "downsizing" trend of a modest evolution of the most 
massive galaxies accompanied by a build-up of small galaxies 
at 2 < 1. We propose that a solution of this discrepancy re- 
quires at least two mechanisms, because the evolution of the 
bright end is driven by inevitable galaxy mergers, while star 
formation in less massive galaxies is sensitive to feedback. 
The solution proposed by Monaco et al. (2006) for slowing 
the evolution of massive galaxies requires an implementa- 
tion of scattering of stars to the diffuse component that is 
strongly dependent on DM halo mass and/or redshift. On 
the other hand, a solution of the overabundance of 10^^ Mq 
galaxies at 2 ~ 1, common to most galaxy formation mod- 
els (Fontana et al. 2006), calls for some unknown source of 
feedback. 
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